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Transverse momentum dependent parton distributions (TMDs) appear in many scatter¬ 
ing processes at high energy, from the semi-inclusive DIS experiments at a few GeV to the 
Higgs transverse momentum distribution at the LHC. Predictions for TMD observables 
crucially depend on TMD factorization, which in turn determines the TMD evolution 
of the observables with energy. In this contribution to SPIN2014 TMD factorization is 
outlined, including a discussion of the treatment of the nonperturbative region, followed 
by a summary of results on TMD evolution, mostly applied to azimuthal asymmetries. 
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1. TMD factorization 

Many angular asymmetries in the transverse momentum distribution of produced 
hadrons in semi-inclusive deep inelastic scattering (SIDIS), ep —>■ e' hX, have been 
measured by the HERMES, COMPASS, and JLab experiments. Evolution is needed 
to compare those results obtained at different energies. The evolution is dictated by 
the appropriate factorization. As SIDIS is sensitive to the transverse momentum of 
quarks through a measurement of Ph±, the observed transverse momentum of the 
produced hadron, it is naturally described within the framework of TMD factoriza¬ 
tion. Several forms of TMD factorization have been put forward in the literature 
for a number of processes [US], which besides SIDIS includes the Drell-Yan (DY) 
process (lepton pair production in hadron-hadron collisions), back-to-back hadron 
production in electron-positron annihilation (e+e" —>■ /11/12A), and Higgs produc¬ 
tion. The main differences among the various approaches concern the treatment of 
spurious rapidity or lightcone divergences, in order to make each factor well-defined, 
and the redistribution of contributions to avoid the appearance of large logarithmic 
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corrections. For brief summaries and comparisons cf. OH] ■ Schematically the TMD 
factorization is of the form |3]: 

da = H X convolution oi AB + high-gr correction (Y) + power-suppressed (1) 


Here A and B are TMD parton distribution or fragmentation functions and H is the 
partonic hard scattering factor. A soft factor has been absorbed into A and B [S] 
|4] . The convolution in terms of A and B is best deconvoluted by Fourier transform. 
More specifically, for SIDIS the differential cross section is given by: 


da 


= j dHe-^'' ’tTW{h,Q-x,y,z) + 0 {Ql/Q^). 


dxdydzdfjid?'Qj, J \ i-xj y ■ j ( ) 

The correction is relevant at large Qt (= Phx!z) and is commonly referred to as 
the F-term. For unpolarized hadrons and quarks of flavor a, W consists of 3 factors: 


W{b,Q\x,y,z) = ^ f^{x,b;CF,y)DUz,b-XD,y)H {y,Q;y). 

a 


( 3 ) 


The Fourier transforms fi and Di of the unpolarized TMD distribution and frag¬ 
mentation functions, are functions of the momentum fraction x or z, transverse coor¬ 
dinate b, rapidity variable and renormalization scale y. Here Cf = 
and Qd = where ys is an arbitrary rapidity that drops out of the 

final answer and CfCf> ~ Q'^i with Q the hard scale. The operator definition of the 
TMDs involves a gauge link or Wilson line U, which arises from summation of all in¬ 
sertions of gluons with longitudinal polarization that are not power suppressed. The 
path of the Wilson lines depends on whether the color flow in the process is incom¬ 
ing or outgoing [ms]. This does not automatically imply that observables depend 
on this path, but it does in certain cases, for example, the Sivers asymmetries |121 
na, where the transverse momentum dependence is correlated with the proton spin 
direction m- The more hadrons are observed in a process, the more complicated 
the color flow, leading to more complicated expressions [mill] or sometimes even 
factorization breaking [T9]El]. In addition, the gauge links in SIDIS and DY have 
lightlike pieces which lead to spurious lightcone divergences. As a regularization, the 
path can be taken off the lightcone, specified by some finite rapidity. The variation 
in this rapidity determines the change of the TMD with C- Also, the regularization 
allows for calculation of the Sivers and Boer-Mulders effects on the lattice |22) . 

Choosing the renormalization scale y = Q avoids large logarithms in the hard 
scattering part i?, but generates them in the TMDs. For this reason one usually 
evolves the TMDs to the scale y^ = Ci/b = jh (Ci ~ 1.123) [3]. This can be 
done using the Collins-Soper and renormalization group equations: 

d\nf{x,b;(,y) ~ dln/(x, 6; C, m) / i \ i ■Z\ (a\ 

- — 7 =- = K{b-y), --- =-fF{g{y);C/y ), (4) 

a In V C a In ^ 

with dK/d\ny = -^K^giy)) & iF{giy);C/= iFig{y)A) - ^7K{g{y))MC/A^)- 
Using these equations one can evolve the TMDs to the scale yb'. 

/r(x, Cf, y) b^-,QD,y) = b^; yl yt) Dliz,b^; yl yb), (5) 
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where the Sudakov factor for Cf = C-D = M = Q is given by 

S{b,Q) =K{b,fib)-[ ^[7F(g(Ai);l)-iln 7ic(5(Ai))]- (6) 

\ Mh / Jfil M ^ \ M / 

The perturbative expression for the Sudakov factor can be used whenever the re¬ 
striction b^ <C 1/A^ is justified (e.g. at very large Q^). If also contributions at larger 
b are important, e.g. at moderate Q and small Qt-, then one needs to include a non- 
perturbative Sudakov factor Snp, for instance as follows: W{b) = IT( 6 *) 
where 6 * = < Vax- For b^ax = 1-5 GeV“\ Q!s(C'i/ 6 max) ~ 0.6, 

such that kF( 6 ,) can be calculated within perturbation theory. In general the non- 
perturbative Sudakov factor is Q dependent and of the form [53]: SNp{b, Q) = 
/Ql)gi{b) + gA{xA,b) -b gs(xs, &), where Qo = l/&max and gi/A/B need to 
be fitted to data. Until recently Snp was typically chosen to be Gaussian, but it 
appears hard to find one universal Gaussian Snp that describes both SIDIS and 
DY /Z production data |5S] . Different b dependences are considered in jS] [HI [55] . 

The TMDs at initial /ii, Q and final /i/, C/ can be related by an evolutor i?, i.e. 
fix, b; Cf, g.f) = R{b; Q,Hi, Q, Hf)fix, b;Ci,g^), with 

i?( 6 ;C.,/i„C/,w) = exp|£'^7F(a.(/l),ln|^)| (7) 

In [27] resummation of logarithms in the perturbative expression for D[b,g) = 
— ^K{b,g,) is performed to NNLL order. To this order the resummed expression 
shows convergence for 6 < 6 x/ 2 , where to leading order bx = ^ exp 

m- The resummed evolutor R vanishes well before b ^ bx/2 ii ^ gi and may 
thus reduce the impact of the nonperturbative b region. Using the 5* method this 
approach favors fomax 1-5 GeV~^ ^7\. Similar resummations in the perturbative 
expansion of the TMDs are performed in [28] . It shows that at low scales the TMDs 
are very small at large b where asifJ-b) is very large. Furthermore, in [28] it is 
suggested that the sensitivity to the Landau pole is minimized by using as initial 
scale Qo -I- qp rather than /if,. Correspondingly a nonperturbative factor with a new 

form is considered: (l + A 2 &^) (Q^/Qq) ^ ^ ■ High Q data (DY/Z) need only 

Ai and A 2 . Low Q data (SIDIS) require modification by including A 3 . 

2. TMD evolution 

Under TMD evolution the shape of the transverse momentum distributions changes. 
Typically TMDs become broader and decrease in magnitude with increasing energy 
[2911,33] . If one starts out with an approximately Gaussian distribution at low scales, 
then a power law tail develops under TMD evolution, see e.g. [31] . This is in contrast 
to a DGLAP evolution of f{x, kp; g) oc /(cc; g) that has sometimes been considered. 
In [23] it was shown that in the limited range of Q from 1.5 to 4.5 GeV, where 
Snp dominates the evolution, such DGLAP evolution hardly modifies the TMDs, 
whereas TMD evolution reduces them by about a factor of 2. 
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TMD evolution has been studied for various azimuthal asymmetries, which gen¬ 
erally decrease as the energy increases. The Sivers asymmetry in SIDIS and DY has 
been studied in [551(5711^ : the Collins effect in e“''e“ annihilation and SIDIS in [551 
IMIISH] ; and the Sivers effect in J/^/> production in [551 SO]- The main differences 
among the approaches are in the treatment of the nonperturbative Sudakov factor 
and the treatment of leading logarithms, i.e. the level of perturbative accuracy. 

First we discuss the TMD evolution of the Sivers asymmetry. The HERMES data 
((Q^) ^ 2.4 GeV^) lie mostly above the COMPASS data ((Q^) ~ 3.8 GeV^) [41]. As 
can be seen in the study of [S^, TMD evolution from the HERMES to COMPASS 
energy scale seems to work well. This result is obtained with some approximations 
that should be applicable at small Q: 1) the Y term is dropped (or equivalently the 
perturbative tail of the TMDs); 2) evolution from a fixed starting Qq rather than /i^; 
3) Gaussian TMDs at the starting scale Qo are adopted. It has been observed in [36] 
that under these approximations plus the assumption that the TMDs as functions of 
are slowly varying functions of b in the dominant b region, the Q dependence of the 
Sivers asymmetry just resides in an overall factor: oc A(QttQ)- Using 

these approximations the peak of the Sivers asymmetry decreases as and 

the peak of the asymmetry shifts slowly towards higher Qt [36] . In [55] it was found 
that the asymmetry integrated over the measured x, z, Ph± range falls off faster than 
1/Q but slower than l/Q^. Testing these features needs a large Q range, requiring 
a high-energy Electron-Ion Collider (EIC). At low (up to ~ 20 GeV^), the 
evolution is dominated by S'tvp [34]. Precise low data can help to determine the 
form and size of Snp, which is responsible for the ±0.1 in CLAS12 is 

projected to have very precise data between 1 and 7 GeV^ (see page 32 of [42]). 

Next we discuss the TMD evolution of Collins asymmetries. The Collins effect 
is described by a TMD fragmentation function [45], giving rise to a sin((()?i ± (j)s) 
asymmetry in SIDIS, in combination with the transversity TMD. Unlike the Sivers 
asymmetry, for the Collins asymmetry no clear need for TMD evolution from HER¬ 
MES to COMPASS (2010) data is apparent. This also needs to be investigated 
further using future data from JLab 12 and possibly EIC. The Collins fragmen¬ 
tation function can be measured independently through the double Collins effect 
(DCE) cos 20 asymmetry in e+e“ —>■ hi /12 X [33], which has been clearly observed 
by BELLE [351I3S], BaBar [37] and BESHI [35]. Under similar assumptions as for 
the Sivers asymmetry, also the DCE asymmetry (and its double ratio for unlike sign 
over like sign hadron pairs) is proportional to an overall factor {Qt) (cf. [49] 1. 
It shows a considerable Sudakov suppression 1/Q [351 |30j |35], which is in rough 
agreement with the results in [251135]. The 1/Q behavior should modify the Collins 
effect based transversity extraction, when full TMD evolution is implemented. Due 
to the lower Q of the BES data, here one does have to worry about 1/Q^ correc¬ 
tions (analogue of the Cahn effect) [501 EHj which can be bounded by studying 
simultaneously the 1/Q cos 0 asymmetry as explained in [30] . 

Finally, we turn to the Higgs transverse momentum distribution, which is also 
a TMD factorizing process. The hard scale Q = Mp is fixed in this case, but TMD 
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evolution matters nevertheless, as the gluon TMDs are probed over a whole range 
of scales fib- The Higgs transverse momentum distribution is sensitive to the linear 
polarization of gluons inside the unpolarized protons [521155] . It requires nonzero 
gluon transverse momentum, but unlike the Sivers and Collins TMDs, it is even 
in kr- Numerical studies of the effects of linear gluon polarization on the Higgs 
transverse momentum distribution vary from permille level [56] to several percent 
m, but with quite some uncertainty from the nonperturbative large-6 region and 
to a lesser extent from the perturbative very small 6 region (6 <C 1/Q). At the 
Higgs mass scale the uncertainty from the latter region is estimated to be less than 
15% by adopting different regularizations (as in Ezl), like the standard one of [58] , 
The treatment of the very small 6 behavior of W becomes increasingly relevant for 
smaller Q values (scalar quarkonium production), and is connected to the treatment 
of the T-term, in order to reproduce the correct integrated cross section, which itself 
should not be affected by the linear gluon polarization. This requires further study. 


3. Summary 

There have been many significant developments on TMD factorization and evolu¬ 
tion recently: new TMD factorization expressions without explicit soft factor and 
with each factor well-dehned have been obtained for several important processes; 
additional resummations have been performed; and, there has been progress towards 
describing SIDIS, DY, and Z production data by a universal non-perturbative func¬ 
tion. TMD evolution has been studied (at varying levels of accuracy) for Sivers and 
(single and double) Collins effect asymmetries and for Higgs production including 
the effects of linear gluon polarization. Future data from JLab 12 and BES and per¬ 
haps a high-energy EIC can help to map out the Q dependence of Sivers and Collins 
asymmetries in greater detail. Future data from LHC on Higgs (and heavy quarko¬ 
nium) production and a high-energy EIC could do the same for gluon dominated 
TMD processes. TMD (non-)factorization at next-to-leading twist remains entirely 
unexplored, but the dependence of azimuthal asymmetries at twist-3 (e.g. 
will be measured in detail at CLAS12, posing a remaining theory challenge. 
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